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Conditions for the existence of strictly stationary multivariate GARCH processes in the 
so-called BEKK parametrisation, which is the most general form of multivariate GARCH 
processes typically used in applications, and for their geometric ergodicity are obtained. 
The conditions are that the driving noise is absolutely continuous with respect to the 
Lebesgue measure and zero is in the interior of its support and that a certain matrix built 
from the GARCH coefficients has spectral radius smaller than one. 

To establish the results semi-polynomial Markov chains are defined and analysed using 
algebraic geometry. 
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1. Introduction 



Generalised autoregressive conditionally heteroskedastic (GARCH) processes (originally introduced 
by (6j [16]]) are heavily used in various areas of applications for the modelling of heteroskedastic 
time series data. Very often one has to model several interrelated time series with an appropriate 
multidimensional model. Since for multivariate GARCH processes the latent volatility process needs 
to take values in the positive semi-definite matrices, as it has to correspond to a covariance matrix at 
each point in time, the multivariate GARCH models are typically considerably more involved than 
the univariate one. For an overview over the various multivariate GARCH models existing and their 
applications we refer to [3TJ . As always in time series modelling, existence and uniqueness of 
stationary solutions, as well as convergence to the stationary solution is of high importance. While 
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this is not so hard a question for some multivariate GARCH specifications where - as in the Constant- 
Conditional-Correlation (CCC) Model or its extensions - the variances are modelled by univariate 
GARCH models, we address this question here for the very general Baba, Engle, Kraft and Kroner 
(BEKK) GARCH model introduced in 1 17] where all (co)variances influence each other in the time 
dynamics. The BEKK model is almost the most general multivariate GARCH model existing. Only the 
vec model also introduced in ifTTl is more general, but all vec models not representable in the BEKK 
parametrisation are somewhat degenerate (see 11321 "). Moreover, the restrictions on the parameters 
necessary to ensure a proper GARCH model are more or less not practicably formulatable in the vec 
model. Hence, the BEKK model is the most general one normally used. 

To prove our results we employ Markov chain theory combined with algebraic geometry to obtain 
the proper state spaces and properties like irreducibility. In particular, we analyse stationarity and 
ergodicity for a general class of Markov chains which we call "semi-polynomial", because they gen- 
eralise the polynomial Markov chains of B71 . and then apply the general results to the special case of 
multivariate GARCH processes. Our approach in the present paper is most similar to the PhD thesis 
of the first author [8 1, which, although it has never been published in an accessible way (except for a 
summary of the main results without proofs in [9]), has been relied upon in essential ways (see e.g. 
ifTOl ). Unfortunately, the statements and proofs in that thesis contain some problematic issues and, 
hence, details of the proofs and statements given in the present paper deviate in essential ways. One 
difference, for instance, is the use of (weak) Feller chains, another regards the proper state spaces. 

The remainder of this paper is structured as follows. Below we briefly summarise some general 
notation. In Section [2] we give a detailed definition of BEKK GARCH models, their vec and vech 
parametrisations, state our main result on the stationarity and ergodicity of multivariate GARCH pro- 
cesses and discuss its implications. Thereafter, we define and analyse semi-polynomial Markov chains 
in Section[3] Finally, we proof our main result for multivariate GARCH processes in SectionlU A brief 
summary of some notions of algebraic geometry necessary to understand the statements of our main 
results on GARCH processes is given in the appendix. There we have also collected some results from 
the theory of Markov chains which we are going to use in the proof of Theorem 13. 121 

Notation 

For the natural numbers excluding zero we write N*. 

We denote the set of real n x d matrices by M„ X ^(R), the vector space of real d x d matrices by 
Mrf(IR), the linear subspace of symmetric matrices by Sj, the positive semi-definite cone by and 
the (strictly) positive definite matrices by For a positive definite and positive semi-definite matrix 
A € we also write A > and A > 0, respectively. The transpose of a matrix A G M„ xt /(IR) will be 
denoted by A' and the d xd identity matrix by 1^. If the dimension is obvious from the context, we 
sometimes neglect the subscript. 

Every matrix A G M„ Xi /(IR) can be considered as a vector in R" rf using the bijective vec transform- 
ation which stacks the columns of a matrix below one another beginning with the leftmost one. In 
the case of symmetric matrices, one often uses the vech transformation which maps §y bijectively to 
jgd(d+i)/2 ky stacking the lower triangular portion of a matrix. For instance, the matrix 



is mapped to the vector (a,b,c,b,d,e,c,e,f)' <E M 9 by the vec operator and to (a,b,c,d,e,f)' 6 M 6 
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by the vech operator. Finally, we denote for two matrices A G M nX d(M) and B G M rxm (M) the tensor 
(Kronecker) product by A 5. 

For the relevant background on Markov chains and mixing we refer to any of the standard refer- 
ences, for instance, |[T3ll26Tl . 



2. Stationarity and Geometric Ergodicity of BEKK Multivariate 
GARCH Models 

When one moves from a single-dimensional to a J-dimensional GARCH process, the univariate vari- 
ance process becomes ad xd covariance matrix process E. In the so-called vec parametrisation (see 
ifTTlO the general multivariate GARCH(/?,g) model with p,q G N is given by 

X n = £ n £ n , (2.1) 
vec(I, ? ) = vec(C) + vec(X n _^_ ( .) + £ Bj vec(E„_ y ) (2.2) 

i=l 7=1 

for n G N* where N* i s an 

Revalued i.i.d. sequence and Z, 1 / 2 denotes the unique positive semi- 
definite square root of £„ (i.e. the unique element l}J 2 of such that L n = l}J 2 l}J 2 ). To ensure the 
positive semi-definiteness of the covariance matrix process £ the initial values £oj • ■ • ^l-p an d C have 
to be positive semi-definite and A\, . . . ,A q as well as B\, . . . ,B p need to be d 2 x d 2 matrices mapping 
the vectorised positive semi-definite matrices into themselves. The initial values Xq,X-i, . . . may 
be arbitrary elements of W 1 . 

The restriction on the linear operators A, and Bj necessary to ensure positive semi-definiteness 
gave rise to the so-called BEKK model (see again [17]) which automatically ensures positive semi- 
definiteness: 

X n = £ ;J £ ;1 , (2.3) 

£ » = C + I I hkXn-iK^ + £ £ Bj^n-jB'j.r, (2.4) 
i'=lfe=l 7=1 r=l 

where A,-^ and Bj r are now arbitrary elements of Mj(R). 

The BEKK model is equivalent to the vec model with A,- = Y!l=\ ^i,k ®M,ki i = !>■•■)?» an d Bj = 
L'/ =1 B j.r^B j_ r , j = I,..., p. More details of the relations between vec and BEKK GARCH models 
are given in Il32l . 

If we take the symmetry of the matrices £„ into account, we can write the vec (and thus the BEKK) 
model also in the vech representation: 

X n = £ n £ n , (2.5) 

q p 
vech(£„) = vech(C) + £ A t vech(X, 1 _ ! -X^,.) + £ Bj vech(L„_ y -) , (2.6) 

i=l ;=1 

where A, = HjAiKj, Bj = HjBjKj, i = I,... ,q, j = 1, . . . ,p and the matrices and y? £ / are the 
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unique H d ,K d £ M jy+i) yd? (R) such that 

vech(D) = H d vec(D), vec(D) = A* vech(D) and H d K f d = I d(d+1)/2 

for every D € whose existence and uniqueness follows immediately, since both vec and vech are 
linear operators. 

Example 2.1. Let us give a simple example for the GARCH equations when d = 2. Consider the 
following BEKK model with p = q = l\ = s\ = \: 



C + 



a c 
b d 



X n — j 



a b 
c d 



+ 



e 8 
f h 



77 — 1 



e f 
h 



(2.7) 



where a,b,c,d,e,f ,g,h are arbitrary real numbers. 

From (12.71 ) one immediately derives the corresponding vec model: 



vec(Z„) = vec(C) + 



/ a 2 ac ca c 2 \ 

ab ad cb cd 

ba be da dc 

\b 2 bd db d 2 j 



vec(X„_ 1 ^_ 1 ) + 



Since the matrices H2 and K2 are given by 
H 2 



1 ON 
10 
,0 1, 



and K2 



(e 2 


eg 


ge 


§2 \ 


ef 


eh 


gf 


gh 


fe 


fg 


he 


hg 


\f 


fh 


hf 


h 2 ) 


1 









1 


1 


I 















vec(E„_!). (2.8) 



the associated vech parametrisation of (I2.71 i and (I2.81 l is 



lac 



leg 



vech(I„) = vech(C) + ( ab ad + bc cd vechpfn-iX,'^) + ef eh + fg 

'■ Ibd d 2 l \f 2 Ifh 



vech(I„_i). 



Definition 2.2. Let (X n ) nE fq* be a GARCH(p,q) process in the BEKK representation satisfying 

(i) C and the initial values Lqj • •• j^l-p are positive definite, 

(ii) (£,i)„ e N* is an ^-valued i.i.d. sequence with distribution T, £\ € L 2 , E[ei] =0 awt/E^eU =I d , 
(Hi) (£ n )nef>r* is independent of = a(Xi^ q ,X 2 - q , . . . ,X ,Ei_ p ,E 2 -/ 

Then (X„) neN * is called standard GARCH(p,<7) process. 



E[ei£j] = I d is a simple normalisation making the volatility process identifiable and, hence, not 
really a restriction. 

1 12 

Remark 2.3. In X n = £„ e n one could also choose another transformation G(L n ) of the conditional 
covariance matrix £„ such that G(r„)'G(£„) = L n instead of taking the square root. Boussama [8] 
takes G(L n ) as a lower triangular matrix resulting from the Cholesky decomposition. However, the 
exact form of this transformation does not have any impact on our results concerning stationarity and 
ergodicity provided that G is an appropriate "smooth" transformation such that the GARCH process 
fits into the setting of semi -polynomial Markov chains studied in Section [3] 
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Using the vech representation we embed a standard GARCH(p,q) process into a Markov chain. 
Setting := (vech(C)',0,0,. . . ,0)' € (lH d+1 )/ 2 ) p x (R d ) q snd defining the process (F„)„ eN * in 

( R d(d+i)/2y x ( E y by Yn := (vech^y^ech^O'^.^vech^+jy,^,^!,...^^ 

one easily sees that 

( Z Ai vech(X„_ ! -X^ i .) + £ Bj vech(E n _ 7 )^ 



-q+\ 



¥„=<<? + 



i=l 



7=1 
vech(S„_i) 



vech(r„_ p+ i) 
X, 
X„_i 



/ 



From this one obtains immediately a regular (in the sense of Definition IA.7I) map 

<p : (Yr^+i)/ 2 Y x (V)'') xR d ^ faW+WY x (r<*) 4 such that Y n = (p(Y n ^,X n ). 

Next we need to define the set W which the stationary standard GARCH(p,q) process takes its values 
in, as is to be seen. Set 

(B\ B2 ... B p _i B p \ 



B : '.o o) GM P^i + ^ M ) withB: 



/ ... 

/ ■■. 



\o ... 



EM 



p 2 



0/ 



and consider the point T satisfying 



T = tf + BT. 



(2.9) 



Existence and uniqueness of this point under the conditions of Theorem I2.4l is shown in Section 
Furthermore, we define cp n for n € N* recursively by q> 1 := <p and 



■bd\n+l 



-l)/2 



p n+1 (y,xi, . . . ,x n ,% n+ i) = <p(9"(y,xi, . . . ,x„),%+i) 



and set 



with S denoting the closure of a set S in the Zariski topology (see appendix). Since L n is always 
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positive definite, we define 

U := vech(S++) x . . . x vech(S+ + ) x R d x . . . x R d 



and consider fW as the natural state space for (Y n ) ne ^*. Finally, HU) denotes the Borel 

a-algebra over W n U inherited from the standard Borel a-algebra on (W l ( d+1 )/ 2 y x (R rf ) 9 , i.e. 
£$(W fl t/) is related to the usual Euclidean topology. 

Theorem 2.4. Lef (X n ) ne fq* be a standard GARCH(p,q) process. 
(0 If 

(HI) the distribution F of E\ is absolutely continuous with respect to the Lebesgue measure on 
R d , 

(H2) the point zero is in the interior ofE: = supp(r) and 
(H3) the spectral radius ofjfj =1 Aj + Y?j=\ Bj lJ i ess than 1, 

then the Markov chain (F n )neN* is positive Harris recurrent and geometrically ergodic on the 
state space (W <~)U ,3§(W DU)). 

The strictly stationary solution (X n ) ne z of the standard GARCH(p,q) model associated with 
(Y n )neZ is unique and geometrically /3 - mixing. Furthermore, X n £ L 2 and L n G L 1 for all 
n £ Z. The strictly stationary solution is, hence, also weakly stationary. 

(ii) If there exists a weakly stationary solution for the standard GARCH(p,q) model, then the spec- 
tral radius of the matrix J^LjAj + Bj lJ ^ ess than 1. 

Remark 2.5. (i) If the initial values Eq, ... are in W and positive semi-definite but not pos- 
itive definite, the geometric ergodicity still holds, since one easily sees that one then has Y p G 

wnu. 

(ii) Obviously, our conditions for strict stationarity imply also weak (or second order) stationarity. 
In the univariate case it is well-known that for a driving noise e with finite variance one may 
well choose the GARCH parameters in such a way that a unique stationary solution with infinite 
second moments exists (see EJ|7J|24l, for example). Extending the above result to cover such 
cases seems not possible at the moment, although in principle the main problem is "only" to find 
an appropriate function for the Foster-Lyapunov drift criterion. 

(iii) Theorem [Z4] should cover most of the BEKK multivariate GARCH models used in applications, 
since one usually wants a finite second moment and uses absolutely continuous noises e (e.g. 
multivariate standard normal or standard ? v -distributed noises with v > 2). Hence, (H3) will 
typically be the only condition that needs checking. 

Note, however, that it is possible to weaken the assumptions (HI) and (H2) to the existence of a 
non-trivial absolutely continuous part of the innovation distribution with zero in the interior of 
its support (cf. |[T3l Section 2.4]). More precisely, if 

(Kl) the distribution F of E\ can be decomposed as F = Fq +F\ with To being non-trivial and 
absolutely continuous with respect to the Lebesgue measure on W 1 (cf., for instance, Il30l p. 
174] where the Lebesgue decomposition of any a-finite measure on ^(R ) is explained), 



6 



(K2) the point zero is in the interior of E = supp(To) 
and (H3) hold, then Theorem [23] remains valid. 

We give a brief summary of the changes that have to be done in the following sections if one 
starts with (Kl) and (K2) instead of (HI) and (H2): in the statements of Theorem l3.ll Proposition 
I3.2l and Theorem 13 .3 l one simply has to replace T by Tq. Further the assumptions (Al) and (A2) 
in Section 13.2.1 1 have to be adapted in the same way as (HI) and (H2) have been weakend. A 
crucial point is to replace the right hand side in Proposition 13.91 by the subprobability measure 
F k (T, ®f =i ro). In the same manner we take v in the proof of Proposition 13. 101 equal to the 
subprobability measure F l {T,® l i=l To). The rest of the paper is not affected. 

(iv) If we completely omitted the assumptions (Kl) and (K2) on the innovations, the mixing result 
might be no longer true (cf. |JT]). However, one could try to extend the idea of [15] where the 
existence of a T-weakly dependent strictly stationary solution for a chain with infinite memory 
has been shown under a Lipschitz-type condition, but the expression of this condition in terms of 
the matrices appearing in the BEKK representation seems to be very delicate. In some cases, like 
for instance the upcoming Example 12.71 the contraction condition is easily verified (in this case 
with the Orlicz function <&(x) =x 2 ). Hence, in such cases lfT31 Theorem 3.1] yields the existence 
of a T-weakly dependent strictly stationary solution of the standard GARCH model even if the 
innovation sequence does not possess an absolutely continuous component. 

Finding an appropriate contraction condition in general appears highly non-trivial but may allow 
to show T-dependence without assuming (Kl) and (K2) and without using algebraic geometry. 
This would imply Central Limit Theorems (CLTs) and the validity of bootstrap procedures (see, 
for instance, 0) for the strictly stationary solution. However, for simulation purposes as well 
as CLTs and Strong Laws of Large Numbers when not starting with the stationary distribution, 
geometric ergodicity and the "right" irreducible state space is very important (see again the up- 
coming Example l2.7l ). The latter seem not to be obtainable under T-weak dependence conditions. 
Moreover, note also that T-weak dependence is a weaker notion than strong mixing. 

(v) Strong mixing conditions are, as T-dependence conditions, a way to derive limit theorems. Any- 
way, the CLT under strong mixing, even if the mixing coefficients decline exponentially fast 
(which is the case for all geometrically ergodic Markov chains), needs a stronger order moment 
condition than the second order ones obtained in Theorem l2.4l (see |fl9Tl and also 112011221 ). In Ifl4l 
(see also ll22l Corollary 3]) it has been shown that for a positive Harris recurrent and geomet- 
rically ergodic Markov chain (X t ) te z on a state space S with stationary distribution n, the CLT 
holds for any real-valued function / defined on S which satisfies f s n(dx)f 2 (x)log + \f(x)\ < °°. 
In our case these references together with the obtained results in particular imply the CLT for 
functions / such that J wnu n(dx)f 2 (x)log + \f(x)\ < C ■ J wnu n(dx)V(x) (where V is the func- 
tion specified in the upcoming proof of Theorem I4.91 l. Hence, for any e > 0, the CLT can be 
applied to (1/2 — e)th and (1 — e)th absolute powers of E„ and X n , respectively. Moreover, com- 
bining our conditions on geometric ergodicity with the fourth order moment conditions of lfl~8l 
immediately gives sufficient conditions for the validity of more classical CLTs. 

One might expect that the set W n U spans the space (R d ( fi,+1 )/ 2 ) / ' x (R d y and that, hence, the state 
space of a stationary GARCH process is "non-degenerate". However, this needs not to be true: 

Example 2.6. Consider the following bivariate GARCH(1, 1) model: 

vech(£„) = vech(C) +A vech(X„_ 1 X^_ 1 ) +B vech(E„_i) (2.10) 
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where A and B are two 3x3 matrices such that the spectral radius of A + B is less than 1 and BA = 0. 
Such a GARCH model can be obtained from a BEKK model with l\ = s\ = 1 and B\\A\\ = 0. 

Starting from the initial point T = (vech(£o) ? 5^o)' gi ven by equation (12.9b . we note that vech(Zo) = 
vech(C) +B vech(£o) and Xq = and obtain by iterating (12.101 ) 



Then W is the Zariski closure of the semi-algebraic set /(M 4 ) (see [4, Theorem 2.3.4]) and W has to 
be strictly contained in M 5 since dim/(]R 4 ) < 4 = dimM 4 . 

Note that the problem of degeneracy in this example lies in the non-invertibility of at least one of 
the two matrices A, B. Indeed, it is easy to see that W is of full dimension and no such degeneracy 
as above can occur if A, B (or An, flu in the BEKK formulation) are both invertible and (HI), (H2) 
hold. 

Moreover one has to be very careful when using GARCH models not to use them outside the state 
space W C\U. Typically one would simulate a stationary GARCH process by starting with an arbitrary 
value and letting the process run. The values are only recorded after a burn-in period. The geometric 
ergodicity ensures that after an appropriately long burn-in period the obtained values can be basically 
regarded as coming from the stationary dynamics. To ensure that this approach works, our results 
show that one needs to start in W H U. One choice of the starting values always possible is T which is 
easily calculated from the parameters by solving a system of linear equations. Let us give an example 
where starting values outside W n U indeed lead to a problem. 

Example 2.7. Consider the set-up of Example 12.6 1 with l\ = s\ = 1 and 



with |a| < 1, \b\ < 1 being non-zero real numbers. Obviously (H3) is satisfied. 

It is then easy to see that the component corresponding to the second variance is constant in W, say 
it equals 022- If we start with an initial value Eq with a second variance (£0)22, then one sees easily 
that (r„)22 = C22 + £ >2 (£«-i)22 for all n G N*. Obviously, this equation has a unique fixed point which 
must be equal to 022 and the right hand side corresponds to an injective map. By induction this implies 
(£71)22 7^ 022 for all ?i G N* if the starting value satisfies (£0)22 7^ 022- 

Hence, for such a starting value Y n W for all n € N and thus the distribution of Y n can never 
converge in the total variation sense to the stationary distribution n. This means that we can never 
have geometric ergodicity when allowing such starting values outside W but in U. 



vech(I, ? ) = vech(I ) +A vech(X„_ 1 X^ 1 ). 



Let / be the regular map from M 4 into M 5 given by 
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3. Stationarity and Geometric Ergodicity of Semi-polynomial 
Markov Chains 



In this section we consider a general class of Markov chains and prove criteria for stationarity and 
geometric ergodicity. We will apply the results later on to the special case of multivariate GARCH 
processes, but the general results of this section seem also of interest of their own, since they should 
be applicable to different models as well. 

We consider Markov chains in W of the form X t+ \ = F(X t ,e t ) where (e t ) te fq is an m-dimensional 
i.i.d. sequence and F is an appropriate map as follows. 

Let V C W be an algebraic variety (cf . Definition IA.5t and U an open subset of W and let F : U x 
W -> U be a C 1 - map such that there exist a C 1 - map / : U x W" -> W and a map <p : W x W -> W 
satisfying: 

(Fl) F{z,y) = (p(z,f(z,y)) for all (z,y) eUx W n , 
(F2) (p(VnUxR m ) cvnu, 

(F3) the map (z,y) <p(z,y) is regular in (z,y) (cf. Definition IA.7I) and 

(F4) for all z £ t/, the map /,(■) = /(z, •) is a C 1 - diffeomorphism from M m onto W and the map 
U x R m — )• M m , (z,y) f7{y) is continuous in (z,y) where f~ (•) denotes the inverse map of 

/z(0- 

The case when F is a regular map, i.e. when f z is the identity, has been considered in [27 ] under the 
name "polynomial Markov chains". Similarly to that paper we use extensively algebraic geometry (see 
the appendix for references and the most relevant definitions) and drift criteria to show the stationarity 
and ergodicity of these Markov chains, but the presence of the additional diffeomorphism f z makes 
all proofs considerably more involved. Moreover, we always need to ensure that we stay in U. 

3.1. Properties of the Image Measure 

In a first step we consider how F acts for a fixed first argument on the noise distribution, which will 
lead to y/-irreducibility conditions in the next section. 

In general, the image of W under F z (-) := F(z, •) is a semi-algebraic set in R" with dimension less 
than n (see [4, Theorem 2.3.4]). Thus the Lebesgue measure of this image is often zero. 

Therefore we need to work with Hausdorff measures (see, for example, lfl2l for a detailed intro- 
duction). We suppose that the algebraic variety V is equipped with a regular measure Hv defined on 
(V,33(V)) where &(V) denotes again the Borel a-algebra over V inherited from the usual Euclidean 
topology. Recall that \iy is said to be regular if, for any A £ 38 (y) and any 8 > 0, there exist an 
open set U G 3B(V) and a compact set K E BB{V) such that K C A C U and /j. v (U\K) < 8. In the 
following we assume that this measure [ly is obtained by equipping the regular set MiV) of V (cf. 
Definition IA.6I ) with an appropriate Hausdorff measure which is extended by zero to the singular set 
y(V) = V\M(V). Moreover we henceforth suppose that T is a measure absolutely continuous with 
respect to the Lebesgue measure on W" with density y. For z € V fl U we denote by F z the image 
measure of T under F z in V fl U. Furthermore we define E := supp(r) which is essentially also the 
domain of positivity of the density y. For the notion of smooth points we refer to lT2"7l A 20] or Il28l p. 
42] and note that the definition makes sense for general C 1 - maps. 

Theorem 3.1. Suppose that zo € V DU and F zo (-) has a smooth point in W n . Then T ZQ is absolutely 
continuous with respect to the measure }Xv an d has support F ZQ {E). 
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Proof. First we denote by Y' ZQ the image measure f Zo (T). One obtains immediately by the Density 
Transformation Theorem that r^ is absolutely continuous with respect to the Lebesgue measure on 
R m with density 



Uy)--=Y{f^\y)) 



1 



IdetD/J/^GO)! ' 



(3-1) 



The support of r' Q is given by f Z0 {E). 

There exists xo G R m such that F zo (xq) = (p zo (yo) is a regular point of V and rank(DF Z0 (xo)) = dim V 
where y = f zo (x ). Since F zo (-) = (p Zo (f Zo (-)), we have DF Z0 (x ) = D(p Z0 (y Q ) ■ Df zo (x ). Since f zo (-) 
is a C 1 - diffeomorphism, the linear map Df zo (xo) is invertible. Thus rank(D(p zo (yo)) = dimV. Since 
(p zo (yo) G &(V) and y G ^?(R m ) (R m is a smooth algebraic variety, i.e. &(W") = R m ), the regular 
map (p zo (-) is smooth at vo and hence dominating (in the sense of ETl A23]). Applying ETl Theorem 
3.1] gives the result. □ 

Proposition 3.2. Let z,o G U. Then, for every s > 0, there exists a > such that 

\T ! Z (B) -T' Z0 (B)\ = \f z {T)(B) - f Zfs {T)(B)\ < e 



for all B G S8(W n ) and every z G t/ with \\z -Zo\\ < a. 

Proof. The image measure V z = f z (T) is absolutely continuous with respect to the Lebesgue measure 
on R m with density yi given by equation (13 - 1 b - 

Let £ > and B G £$(W n ). The space of real-valued continuous functions on R m with compact support 
is dense in the L 1 sense in the space of all Lebesgue-integrable functions on R m . Thus, there exists a 
continuous function y : W — > M with compact support K such that 



Hence, 



\r' z (B)-r zo (B)\ 



< 



\y(x)-y{x)\ dx< 



Y z (y)dy- I Y zo (y)dy 

B J B 

7(f-\y)) - ?(f-\y))\ ■ \detDf-\y)\ dy 



(3.2) 



+ 



dy 



■ / f(f-\y))\detDf- l (y)\-f(f- \y))\detDf- \y)\ 
+ j R jY(f^(y))-Y(f- \y))\-\d & tDf- o l (y)\ dy =: h+h+h- 

With (13.21 ) we obtain immediately by substitution I\ < e/3 and It, < e/3. 

y is bounded on R m by sup y and for all r > such that B(zo,r) := {z G R" : ||z — Zo I < r} C £/, the 
set C := { (z,y) G £/ x R m : \\z — Zo\\ <r and /~ (y) G A'} is a compact set in U x R m , since the map 
y/ : U x R m — > U x R m , V'foy) : = (z>/z60) is continuous and C = y(B(z.o,r) x X"). Thus, there is a 
real number b > such that, for all (z,y) G C, | detD/^ 1 (y) | < b. The map y(/ z _1 (y)) ■ | detD/^ 1 (y) | 
is hence bounded on C by Z? • sup y. 

Let Ci be the projection of C on R m and suppose without loss of generality \\z — Zo|| < t. Then, for all 
y ^ Ci, we have y{f z { (y)) = y(/?o 1 60) = which implies 



h= I nf z \y))\teWf-\y)\-%f z -\y))\teWf zo \y)\ 



dy. 
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This integrand is dominated by 2b sup y and converges pointwise to zero if z converges to zo (cf. (F4)). 
Since b and sup y are finite constants and Ci is compact the dominant 2Z?sup y is integrable over Ci. 
Hence, we can apply the Dominated Convergence Theorem and get I2 — > as z — > Zo, i.e. there exists 
< a < r such that 7 2 < e/3 for all z G E/ with ||z - Zo || < «. □ 

Theorem 3.3. Suppose that zq G V fl J7 a«<i F zo (-) /j^s « smooth point in M m . Then 

limMr z (A)>r zo (A) (3.3) 

zevnu 

for every A G 3S{V fit/) (?/je Bore/ a -algebra inherited from the usual Euclidean topology). 

Proof. Let A G &(V DU) and £ > 0. Since <p zo (-) is dominating and FL = / zo (r) is absolutely con- 
tinuous with respect to the Lebesgue measure on W" (cf. proof of Theorem 13- lb . 11271 Theorem 3.2] 
yields a neighbourhood Vq of Zo in V n C/ such that <p z (r' )(A) > <p Zo (r' )(A) — f for all z G Vb which 
is equivalent to 

^(n^-HA)) >/ zo (r)(<p zo 1 (A))-|vzGVo. (3.4) 

Due to Proposition O there exists a > such that / Z (T)(B) > f Z0 (T)(B) - f for all G 38(W) and 
every z G U with ||z — zo || < OC. We choose B = (p~ l (A) and deduce for every z£U, \\z — Zo\\ < OC, 

MT)((p-\A)) > fzo (r)(<p-\A))-~. (3.5) 

With (fill) and 63 we obtain for all z G V with ||z - z 1 1 < a that / t (r) ( q>~ 1 (A) ) > f zo (T) ( p~ 1 (A) ) - 
£. Since <p z (r^) = F z (r) = T,, this is equivalent to 

r z (A) >T Z0 (A)-£ VzGV n{zGM": Hz-zoll <a}. 

This shows (13.31 ) since £ > can be chosen arbitrarily small. □ 

3.2. Stationarity and Ergodicity 

Culminating in Theorem 13.111 we now gradually show Harris recurrence, geometric ergodicity and 
j3 -mixing for semi-polynomial Markov chains. 

3.2.1. Assumptions 

Concerning the sequence (e t ) t eN we make the following additional assumptions for our semi-poly- 
nomial Markov chain: 

(Al) Every e t has distribution Y which is absolutely continuous with respect to the Lebesgue measure 
on W with density y. Let E denote the support of T. 

We define for all k G W,k > 1 the functions F k (z,yi, . . . ,yt) := F(F k ~ 1 (z,y ]_,..., yk-i),yk) where 
ZGU, (y u ... ,y k ) G (R m ) k and set F l = F. 

With this notation we introduce for z G V fl U the orbit 

Sz-= U {F k (z,yi,...,y k ):yi,...,y k eE}. 

ken* 

To prove the desired properties for semi-polynomial Markov chains we assume: 
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(A2) There is a point a G int(is) and a point T £V HU such that, for all z G V n £/, the sequence 
(X, z ) ?e N defined by Xfi = z and Xf = F(Xf_ v a) for t > 1 converges to the point 7\ 

T is called attracting point of the chain (X f ) /e ^. 

We set W : = z 5r the Zariski closure of the orbit Sj- Note that obviously T G W and IV C V. To 
show uniqueness of the strictly stationary solution we need the assumption: 

(A3) Any strictly stationary solution of the Markov chain X t+ \ = F(X t ,e t ) takes its values in the 
algebraic variety W H U. 

Remark 3.4. (i) If (A2) is satisfied, then T is a fixed point of F(-,a), since F is continuous. 

(ii) It is obvious that W is an algebraic set since it is the Zariski closure of St- In fact, it is even 
irreducible (cf. the upcoming Section |3.2.2| |. 

Strictly speaking W n U is not necessarily an algebraic variety, but, as it is the intersection of an 
algebraic variety in W and the set U where our Markovian dynamics are defined, we refer to it as an 
algebraic variety. 

3.2.2. Algebraic Variety of States 

In this subsection we suppose that the assumptions (Al) and (A2) hold. We will show that W, defined 
as above, is indeed an algebraic variety which we will call the Markov chain's algebraic variety of 
states. 

Let (AOfceN* be the sequence of subsets of U defined by D k := F k (T,E k ). Since F(T,a) = T 
(cf. Remark [33](i)) we obtain D k = F k (F(T,a),E k ) = F k+l (T,{a} x E k ) C D k+1 , i.e. the sequence 
(D^keN* is an ascending sequence of subsets of R". 

We set Wk := z F k (T,(R m ) k ). Then we have W k = z (p k {T, (R m ) k ) (defining (p k analogously to F k ) 
since /r(-) is a C 1 - diffeomorphism. 



Proof. To this end consider the map fj : (R m ) k -4 (R m ) k , (yi,... ,}>k) ^ (xi,... ,xt) where x\ 



Lemma 3.5. For all k G N* we have Wk = z Dtc- 

Proof. To this end consider the map /p : (W 

fr(yi), X2 = f F (Ty 1 )(y2),---,x k = f Fk -i {T j u ^ yk l) (y k ). 

(k) 

Due to the properties of / and F (in particular (F4)), it is clear that f T is bijective, continuous and its 
inverse is continuous as well, i.e. ff' is a homeomorphism. 

Assumption (A2) implies that E k contains an open ball of (M" 1 )*. Thus, since f^ is homeomorphic, 
ff'{E k ) contains an open ball of (M. m ) k . 

From H Corollary 3.4.5] we obtain f T k) (E k ) = (R m ) k . This shows 



W k = z F k (T, (W"f) = z (p k (T, (M. m ) k ) = V (t^/tHe*)) ■ 

Since (p k (T,-) is regular (cf. (F3)), q> k (T,-) is continuous with respect to the Zariski topology due 

to Proposition |AU Hence, W k = «p k (r,ffl = z F k (T,E k ) = ^ which proves W k to be the 
Zariski closure of D k . □ 

Lemma 3.6. W k is irreducible for all k G N*. 
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Proof. If we suppose that there is k £ N* such that W% = V\ U V2 where V\ and V% are algebraic sets 
with V\ C and V 2 £ W*, then 

= (w k ) = (Vi)u (y 2 ) 

(*) (**) 

where <?>£(•) = <p fe (r, •). Now (*) and (**) are algebraic sets, because V\ and V 2 are algebraic sets and 
<Py(-) is continuous with respect to the Zariski topology (see proof of Lemma l3.5l >. Since (<p T ) 1 (Vi) C 
(W n ) k for i = 1,2 (otherwise = z <p i (r, (R'")*) C z V t = V t which would be a contradiction to 
Vi C Wyt), this would prove (W) k to be reducible which is a contradiction. □ 

Proposition 3.7. There exists I £ N* smc/j f/ia? W/t = Wifor all k>l and W = Wj. In particular, W is 
an algebraic variety. 

Proof. From [4, Corollary 3.4.5] it follows that if V\ C V% C V3 C ... is an ascending sequence of 
algebraic varieties in R" then there exists / £ N* such that Vk = Vi for all k > I. 

Lemma 1331 and Lemma I3T61 show that (W k ) ke ?$* is an ascending sequence of algebraic varieties and 
so there exists I £ N* such that W k = W/ for all k > I. We then observe that 

St= \J {F k (T, yi ,...,y k ): yi ,...,y k eE}= (J D k . 



ken* s v ' feel 

-F k {T,E k ) 



Since 



^ [j D k Q |J z D k = z W,=W t 

feeN* feeN* 



= U w k =w, 

keN* 



and Wi = z Di C z U /te N *D*, we obtain IV = z S r = W/. □ 

Lemma 3.8. For all k £ N* we have F k (W n [/, (R m )*) C IV n U. Hence, the Markov chain can be 
restricted to the variety of states W n U. 

Proof. With the definition of the subsets D k and W k , respectively, one has for all k £ N* 

(p(D k ,R m ) = F(D k ,R m )=F( F k (T,E k ) ,R m ) C F k+X (T,(R m ) k+l ) C nt/ C?nl/. (3.6) 

cF*(r,(M m )*) 

The continuity of regular maps with respect to the Zariski topology and the regularity of cp yield 



-ED 



F(ffnf/,l m ) = <p(Wn£/,IR m ) C (p(W,M. m ) = (p( z Dt,R m ) C <p( z D,,IR'») = z <p(D,,R" ! ) C W. 

Since we assume F : U x R m ->■ £/, we have F(W n £/,R m ) C W n £/. By induction we have F*(W n 
£/, (R m )*) C W fl C7 for all & £ N*. Hence we can restrict the Markov chain to the variety of states 

wnu. □ 
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Proposition 3.9. For all A G £%(W HU) and allk>l 

liminf^(z,A) >P k (T,A), 

Z-^T 

zewnu 

where P k is the k-step transition probability kernel of the Markov chain (X t ) t€ ^. 

Proof. Since (p T (-) = (p k (T, •) is regular and dominating for all k > I (since z (p k (T, (R m ) k ) =W k = W 
for all k > I, cf. Proposition 13.71) . E71 A 23] implies that <p*(-) has a smooth point. 

Similar to the proof of Theorem 13. II we can now show that the map F k (T,-) has a smooth point in 
(W") k . Let x G (R m )* be the smooth point of <p£(-), i.e. (p k (T,x Q ) G ^(W) and rarajfc (Z><p£(x )) = 

dim W. Then F*(V, f/f^ 1 (x )) = (p k {T,x ) and 

DF k (T, •) ( (/f ) ' (x )) = Dcp k (x Q ) ■ Dff ( (x )) 
where the linear map Df T k \x) is invertible for all x G (R m )* (to this end note that ftp is not only 

(k) 

continuous but also differentiable and that Dfj- is a block matrix with lower triangle structure where 
the blocks on the diagonal are invertible). Hence the matrix on the left hand side has also rank dim W 
and (fx) 1 (xo) is a smooth point of F k (T, •). 

Finally, note that P k (z,A) = F k (z,® k =l T)(A) where F is the distribution of every e t (cf. (Al)) and we 
conclude with Theorem [33] □ 

3.2.3. Harris Recurrence, Ergodicity and p - Mixing 

In this subsection we will prove the promised properties of semi-polynomial Markov chains under a 
Foster-Lyapunov-condition. First we show irreducibility on the algebraic variety of states and aperi- 
odicity. As usual, y denotes a maximal irreducibility measure. 

Proposition 3.10. Suppose that (Al ) and (A2) hold. Then the semi-polynomial Markov chain (X t ) te ?q 
is \\f - irreducible and aperiodic on the state space (W D U,,^(W (~)U)). 
Moreover, the support of \jf has non-empty interior. 

Proof (1) Due to Proposition [3791 we have for all A G &(W DU) 

liminfP'(z,A) > P'(T,A). (3.7) 

zewnu 

We define a probability measure v on the state space (W n U, 3S(W n U)) by v(A) := P ! (T,A),A G 
3S{W PI U). Then, for every A G 3S(W n U) with v(A) / 0, there exists due to (13771 a neighbourhood 
W\ of T in W n U such that 

P'Ml^VzeffL (3.8) 

(2) Let K = {zi, . . . ,z r } Q W D U for some r G N*. We are going to show that there is a q G N* such 
that P q (zi , Wi ) > V/ G { 1 , . . . , r} . To this end, consider for i = 1 , . . . , r the sequences (X t Zl ) te ^ defined 
by 

Xj? = Zi and Xf' = F(X®_ 1 ,a), t > 1 
where a G int(£) as in (A2). 
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Due to assumption (A2) there is q G N* such that X® = F q (zt,a, . . . ,a) G W\ Vi € {1, . . . ,r}. Since 
: W n C7 x (R m ) 9 — > W PI £/ is continuous, there exists for every i G {1, . . . , r} a neighbourhood £/; 
of (zj,a, ...,a)inWnf/x (R m ) q such that 

^(y,yi,.-,y 4 )GWiV(y,y 1 ,...,y,)€yi. 

Then, for all i G { 1 , . . . , r}, Ui contains U- x Ui n where £// and Ui n are suitable neighbourhoods 

of n in W PI U and (a, . .. ,a) in (M m )^, respectively. 

r 

We define U^ a ^ a ^ := P| E/i ^ which is clearly also a neighbourhood of (a,..., a) in (M m ) 9 . 

Then we have F q (z i: yi,...,y q ) G Wf for all i G {l,...,r} and (yi,...,y q ) G £/( a ,..., a ). Since t/^...,,,) 
contains itself U a x ... xU a where C/ a is an appropriate neighbourhood of a in M m , we deduce for all 
ie{l,. ..,r}: 

P^Wi) >P((ei,...,e,) G a) ) >P(«i G U a f = T{JJ a ) q . (3.9) 

(3) Let A G ^(W n C/) with v(A) / 0. As in (1), Wi denotes the neighbourhood of T in IV n U such 
that P'(z,A) > v(A)/2 for all z G Wi. Using the Chapman-Kolmogorov equation (cf. E6l Theorem 
3.4.2), we obtain for every / G {1, . . . ,r} 

P<? +/ (z ; -,A) = / P^fc,^)P'(3;,A) > / P?( Zl ,^)p'( y ,A)?^ - / P^rfy) 
v(A) ED r(t/„) 9 

= m./«( a>Wi ) > L -^--v(A). 

Due to assumption (A2) U a contains an open set of E. Thus T{U a ) > 0. This implies that the chain 
(X f ) f6 N is V - irreducible (and thus also \j/ - irreducible due to ll26l Proposition 4.2.2]). 

(4) To show aperiodicity, we suppose the chain to be periodic with period d. Due to |[26l Theorem 
5.4.4] there exist disjoint sets Di,...,D d £ 3§(W n U) such that 

(i) P(z,D (! - modd)+l ) = 1 Vi = 1,. . . ,d andz G D t and (ii) ^((uf =1 D,-) c ) = 0. 

Since \ff >- v (cf. IT261 Proposition 4.2.2]), (uf =1 D,) c is also a v - null set. Obviously there must be 
a set Dj with positive v - measure, let this set be D\ without loss of generality. 
Let x G Di and y G Dj. For := {x, y} we have just shown in step (3) that 

P 9+ '(jc,Di) > and P« +/ (y,Di) > 

for some q G N*. Hence, the integers q + l and g + Z — 1 are divisible by d. Consequently d = 1. 

(5) We have shown in step (3) that {X t ) te m is P l {T, ■) - irreducible. Since F\T, ■) has a smooth point 
in (M. m ) 1 (cf. proof of Proposition 13 .9t Theorem 13. 1 l implies that P l (T, •) = F l (T, <S>- =i r) is absolutely 
continuous with respect to the measure jlw Hence, int(supp P'(r, •)) / and we also obtain that 
int(supp y) ± 0. □ 

We can now state our main result for semi-polynomial Markov chains. Therefore we use the stand- 
ard notation PV(x) := E[V(Xi)|X =x) = E x [V(Xi)]. 

Theorem 3.11. Suppose (Al ) and (A2) are valid. If in addition the Foster-Lyapunov -condition holds, 
i.e. there exist a small set C G 38 (W DU), positive constants a < 1, b < °° and a function V > 1 such 
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that 



PV(x) <a-V(x) + b-l c (x)VxeWnu, 



(FL) 



then the semi-polynomial Markov chain (X f ) fe pj is positive Harris recurrent and geometrically er- 
godic on the algebraic variety of states W D U. Furthermore, the strictly stationary process (X t ) te i is 
geometrically j3 - mixing and 7l(V) := K[V(X t )] < oo. 



Proof. Due to Proposition 13. 101 and the assumptions (Al) and (A2), (X t ) te ® is \j/ - irreducible and 
aperiodic on (W n U,^(W H t/)). We conclude by using the following Theorem 13. 121 □ 

Theorem 3.12. Let (X t ) te fq be ay - irreducible Markov chain on a state space (S,£$(S)) with trans- 
ition probability kernel P. If the chain is aperiodic and the Foster-Lyapunov-condition holds, i.e. there 
exist a small set C G £8(S), positive constants a < 1, b < oo and a function V > 1 such that 

PV(x) <a-V(x)+b-\ c (x)Vx£S, 

then (X t ) te fq is positive Harris recurrent, geometrically ergodic and the strictly stationary process 
(X t ) t £Z is geometrically /3 - mixing. Furthermore, lt(V) < °°. 

Proof. Since the Foster-Lyapunov-condition holds, the non-negative functions V' := V — 1, / := 1 — CC 
and s := blc satisfy the assumption of Theorem IB. II Hence we obtain for the first return time to C, 
denoted by Tc, 



1 



1 - a 



E v 



k=Q 



< 



1 



1 - a 



7(jc)-1+E, 



I s(X k ) 
k=Q 



< o°Vx G S 



=bl c (x) 



and thus obviously L(x,C) = f x (^c < °°) = 1 for all x G S. Since every small set is also petite (cf. [25 ] 
or 112610. Proposition IB .21 yields Harris recurrence of (X t ) te fq. Again the Foster-Lyapunov-condition 
shows that V' := (1 — a)~ l V, f := V and b' := (1 — a)~ l b satisfy (ii) of Theorem IB . 3 1 which implies 
that (X t ) te ® is positive and 7t(V) < °°. It is once more the same condition that yields directly geometric 
ergodicity by virtue of Theorem IB .41 Finally, combining iTTTl Proposition 1 (1)] with (IB.ll ) and 7t(V) < 
oo, we deduce that the strictly stationary process (X t ) te z is geometrically /3 -mixing. □ 



Theorem 3.13. Suppose the setting ofTheorem \3.11\ and assume in addition that (A3) holds. Then the 
strictly stationary process is unique. 

Proof. If there is another P - invariant probability measure %' , then supp(7r') C W P\U due to (A3). 
Since the chain (X t ) te ^ is recurrent on W n U (Theorem 13. UK it has at most one P - invariant prob- 
ability measure on (W D U,88(W D V)) (cf. ll26l Theorem 10.4.4]). Therefore the strictly stationary 
solution is unique. □ 

Remark 3.14. Note that the whole theory from algebraic geometry has only been used to prove irre- 
ducibility and aperiodicity on the algebraic variety of states W n U. The results thereafter (Theorems 
13. 12l and l3. 13b are consequences of the theory of Markov chains. 



4. Proof of the Main Theorem 2.4 



In this section we gradually prove our main result for multivariate GARCH processes, Theorem 
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4.1. GARCH Processes as Semi-Polynomial Markov Processes 

First we show that the autoregressive representation of standard GARCH processes involving the 
function (p leads to a semi-polynomial Markov chain. 



++ 



Lemma 4.1. The mapping G : St — > E>^ which maps £ to Z 1 / 2 is aC l - diffeomorphism on S 

Proof. We define f d : §++ -> §++, X ^ X X. Since every X € S++ has full rank, / d is well-defined. 
We will show that fd is bijective and that, for every X € §>,} + , the differential df d (X) is a linear 
homeomorphism. 

By ETl Theorem 7.2.6], a positive definite matrix has a unique positive definite square root. Hence 
fd is bijective. Let X G The differential of at the point X is given by 



Vtf € S d df d (X)H = ^f, (X + tH) 



= HX+XH. 

f=0 



Our aim is to show that H = whenever dfd{X)H = 0. In fact, this is a simple consequence of 
Theorem 1] where the solutions X of the general matrix quadratic equation = A + BX +XB' —XCX 
for fixed A,B,C € Md (R) have been analysed. □ 

1 12 

Since X„ = L n s„, we thus obtain that 

f:UxR d ^R d , (Y n _ l ,s n )^G(L n )s n =X n 

is a C 1 - map from U x R d into M rf where U is the open set in (R^+i)/2) p x defined by 

U := vech(S+ + ) x . . . x vech(S+ + ) x R d x . . . x R c [ . 

p * 

Due to the assumption that C and the initial values Lq, . . . ,£i_p are positive definite, every £„ and 
Z, 1 / 2 is also positive definite and thus G(L n ) = l}J 2 is always an invertible matrix. Hence, for every 
Y G U, the map fy (•) = /(Y, •) is linear bijective from R d onto R d , i.e. /y(-) is a C 1 - diffeomorphism. 
Moreover the map U x R d — > R d , (7,£) h-> fy (e) is continuous in (F,e) where fy (•) denotes the 
inverse of /y(-). 

Altogether we are thus in our setting of semi-polynomial Markov chains. For the Markovian rep- 
resentation Y n of a standard GARCH(/?,g) process X„ we have 

Y n =F(Y n ^,e n ) := <p(Y„_i,/y n _, (e„)) (4.1) 



where T 7 is a C 1 - map from U xR d into £/. Moreover, it is obvious that (I2.5I ). (I2.61 l have a stationary 
solution if and only if (14.11 ) has one. 

4.2. Some Results from Linear Algebra 

In this section we will show some results from linear algebra which will be necessary to establish the 
Foster-Lyapunov-condition for multivariate GARCH models. 



17 



Let G N* and (F/)i</<„ be elements of M d (R). We set 

£ :M d (R) ->■ Afd(M), £(M) := £f'MF/. 

This map is obviously linear. We can consider t, a linear map from R rf2 into M rf2 using the vec operator 
as follows : 

(n \ n 

^F, ® F; J vec(M) = Fvec(M) where F := £f ; - <g>F;. 
(=1 / ;=1 

Note that we have <^ (S^) C § rf , i.e. the symmetric d xd matrices are mapped into themselves by t, . 
We denote by t, the restriction of t, to the linear subspace S d . Using the vech operator, we obtain, for 
all M G S d , 

vech (| (M)) = vech(£ (M)) = # d vec(£ (Af)) = # d F vec(M) = #<*F^ vech(M). 

Since we can identify S</ via the vech operator with ]R rf ( rf+1 )/ 2 , the transformation matrix of | is given 
by 

F:=H d FK' d . 

We obtain the following lemma: 

Lemma 4.2. Le? C € The following statements are equivalent: 
(i) The spectral radius oft, is less than 1. 
(7/J 77ie spectral radius oft, is less than 1. 
( iii) There is E € §j~ + swc/i E = C + £ (E). 

Proof, (i) (ii): Obvious since <^ is a restriction of t,. 

(ii) => (iii): If the spectral radius of t, is less than 1, then the Neumann series £~ = o£" is convergent 
with respect to a suitable operator norm. We define 

CO 

E:=£|"(C). (4.2) 

n=0 

Clearly, £ is symmetric. Moreover, for all M G St, we have t, (M) G 8^ by the definition of t, and £ , 
respectively. By iteration we obtain that | n (M) G for all n G N*. Thus the matrix E-|°(C)=E-C 
is symmetric and positive semi-definite. This implies that E is positive definite. 

Since § and f coincide on § d , we deduce E = E§"(C)=C + q £ =C + €(£). 

«=0 ^n=l ' 

(iii) =>- (i): Suppose that there exists E G such that E = C + t, (E). We denote the complex d xd 
matrices by M d (C) and the conjugate transpose of a vector x G C d by x*. For every P G Af^(C) we 
define 

j|F|| z := sup \x*Px\ 

xeC d , x*T.x=l 

which is a norm on M d (C) since E G + . 
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Then, for all x G C d , \x*Px\ < \\P\\^ (x*Lx). Since the unit sphere {i £ C* : x*Lx = l} is compact, 
there exists, for every P G Af^(C), a vector x p G C rf such that ||P||j; = |xp,Pxp| withxpZxp = 1. Let 

n 

now A be an eigenvalue of % . Then there is an M G M rf (C) , M 7^ such that AM = £ (M) = £ F t MF[ . 

1=1 

For every 1 S C 1 ', we deduce 



|A| • \x*Mx\ 



j^x*FiMFlx < £ I {F?x)*M (Fix) \ < ||M|| E = ||M|| e jc* *. 

i=l i=l i=l \i=l / 



If we choose xm such that \\M\\ Z = \x* m Mxm\ and x* m 1jcm = 1, we obtain |A| < 1 — x* m Cxm < 1 (note 
that ||M|| E 7^ 0). Hence the spectral radius of E, is less than 1. □ 

We consider now the families of matrices (Ai^,Bj r ) and (Aj,Bj) which occur in the BEKK and 
vech representation of a standard GARCH(/?,g) model, respectively. 

Proposition 4.3. The spectral radius of I^LjA, + £jLi^ 'j is less than 1 if and only if there exists 
£ G §j + such that 

1 = C + 1 1 A.v£Aj,t + £ £ tf;,^,.. (4.3) 

i=U=l y=l r=l 

Proof. Define £ :M d (R) — >■ Af^QR) by t, (M) = £ £ A ; - fc MA^ + £ £ B^MB'.. Then the trans- 

i=lfc=l ' y=l r=l 

formation matrix of £ is 

9 k _ p s i q p 

f = £ £ ^ ® 4* + L L ® s i> r = + E s j- 

i=U=l 7=1 r=l i=l 7=1 

Note that the transformation matrix of E, (restriction of £ to the linear subspace is HjFKj = 
Lf=i M + Ey=i ^7- Due to Lemma 14721 the spectral radius of £f_i A, + Yl]=\ Bj is l ess than 1 if and only 
if there is £ G § j + such that @3]> holds. □ 

Remark 4.4. By a simple transposition argument one can equivalently state that the spectral radius 
of A, + Lj=i ^7 i s l ess than 1 if and only if there exists £ G such that 

1 = C + E I ^U 5 ^* +11 S )^Bj,r 
i=l*=l 7=1 r=l 



which we are going to use in the upcoming proof of Theorem [ 

In the following we consider the block matrix B defined as in Section [2] 

Proposition 4.5. (i) If the spectral radius of the matrix Yfj=\ Bj is less than 1, then the one ofB is 
also less than 1. 

(ii) If the spectral radius of the matrix J^LjA, + Y?j=\ Bj is less than 1, then the one ofYfj—iBj is 
also less than 1. 
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Proof, (i) Suppose that the spectral radius of Y?j = \_Bj is l ess than 1. Then there exists due to Lemma 
4.21 a symmetric positive definite matrix £ £ such that 



7=1 r=l 



(4.4) 



Let X be an eigenvalue of B associated with the eigenvector h = (h\, . . . ,h p y £ (M. d( - d+l ^ 2 ) p . Then 
Xh\ = Ylj=\Bjhj an d Xhj = hj-\ for 2 < j < p. Thus /j p ^ (otherwise h would be zero) and 
XPh p = X(XP- l h p ) = Xhi = Yfj =1 Bjhj = Ey=i XP-jBjh p . LetM € S rf such that vech(M) = ^. Then 
X P M = Y?j=\ Lr=i W~ 3 BhrM& j r . We define the norm ||-|| £ on M d (C) as in the proof of Lemma K 
by||P|| £ := sup \x*Px\ , P € M^(C). Then, for all x e C rf , 



|/l| p -|x*Mx| 



£ l^XP-^rMB)^ 

j=\r=\ 



7=1 r=l 



<MfiEElArV*^^)- 

7=1 r=l 

If we assume that there is an eigenvalue A of fi with | A | > 1, then we obtain, taking the vector x such 
that;c*£x = 1 and \x*Mx\ = ||M|| £ and using (14.4b . that 



\M P < 1 1 |Ar^(x*y5 ;> I75^) < \xr i 

;=1 r=l 

= |A| P " 1 [x* (£-C)jc] =|A| p " 1 (l-x*Cx). 



Since C is symmetric positive definite, one has x*Cx > 0. Hence, \X\ P < i.e. |A| < 1 which is 

a contradiction. Thus the spectral radius of B has to be less than 1 . 

(ii) Suppose that the spectral radius of the matrix Ef = i A,- + EjLi Bj is l ess than 1. Then, due to Pro- 
position S3] there exists I G § ++ such that E = C + Ef= j EL [ A i>k TA\ k + lf j= l E*i t Bj^B'j. r ■ We set 
C := C + £f =1 ELi^i^i,*- Now > ^ is symmetric positive definite and I = C + EjLi Erii Bj r LB'j r . 
Using again Proposition 14 . 3 1 we deduce that the spectral radius of Ey=i ^7 * s l ess th an 1- D 

Remark 4.6. The matrix £ in (14.41 ) is the limit of a Neumann series (cf. proof of Lemma l4~2l) .Thus 



vech(£) = ( 7 - L S j) vech (C). 



4.3. Verification of Assumption (A2) 

We will suppose throughout that (HI) holds. 

p 

Proposition 4.7. Suppose that (H2) holds. If the spectral radius of the matrix E ^7 w than 1, 

7=1 

f/jen (A2J /joWs. 
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Proof. Let U be the open set in (]R d ( rf + 1 )/ 2 ) / ' x (R d ) q denned as in Section I47T1 For arbitrary y G U 
we define the sequence (Y^) n€ ^ by Yq =y and Y% =F(7^_ 1 ,0), n> 1. 

We denote by X„ and vech(E^) the associated values of X n and vech(£„). Since, by definition, 

1 /2 

X„ = G(L n )s n = £,/ £„, we obtain that X„ = for all n > 1. Due to (12.6b . can be written, for every 
n > q, as 

vech(£>,) = vech(C) + £ B ; vech(L£_ .)• 
Thus, for all n > q and for all .y G t/, 

Yy = V + BY y n _ x (4.5) 
where B is defined as in Section |2] Due to (14.51) . the assumption (A2) is satisfied with a = if the 
spectral radius of B is less than 1 . This is the case, since the spectral radius of £ Bj is supposed to be 

7=1 

less than 1 (cf. Proposition 14.51 (D). □ 



Hence, for all y G t/, the sequence (F« ) n€ N converges to the unique fixed point T defined by 

T = tf + BT. (4.6) 

Using Lemma [4721 and the fact that the spectral radius of E/=i Bj * s assumed to be less than 1, there is 
£ G Sj^ + (cf. (14741 ) such that 

;=lr=l 

It is then easy to see that T can be written as 

T = ^vech(£) f ,...,vech(£)',0 :: ^0^ eU. (4.7) 

P 

We set <Ti := (vech(C)' ,0, . . . ,0)' G (M^ 1 )/ 2 )''. Then @3 yields 

o = %+Bo (4.8) 
where a := (vech(£) f , . . . , vech(£)') f G (R^ +1 V 2 ) P 

4.4. Verification of Assumption (A3) 

If (H2) is satisfied, then E contains an open set of M. d and we obtain for the algebraic variety of states 
with the same arguments as in Section 13.2.21 that 



W = Z S^= |J F n {T,E n )= |J F n (T,(R d ) n ) = \J (p n (T, (R d ) n ) . 

hgN* neN* neN* 
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Let neW and consider y(n) € <p" (T, (M d ) n ) given by y(n) = <p n (T,x\ , . . . ,x n ) where x\, . . . ,x n £ Mr. 
We define x(n) and o(n) by the coordinates of y(n) as follows: 

x(n) = (vech(x„4) r ,--- 5 v ech(x„_^ +1 x^ +1 )')' 
and a(n) = (vech((T„) f , . . . , vech(<T„_p+i) ? ) . 



That is, y(n) = {o(ny ,x* n , . . . ^ n _ q+l )' ■ Then 

o(n+\) =^i+Ax(n)+Bo(n) 
where ^ and B are defined in Section 1431 and A is given by 



(4.9) 



/Ai A 2 








€M rf(rf + i) 

F 2 x " 2 



\0 ... 0/ 

Iterating (14.91 ) and due to a(0) = a = (vech(£)', . . . , vech(Z) f ) ? (since y(0) = T) we deduce 



n i 

a(n) = £ B'Vi + B"a + £ B'^A^n - i) = a + £ ff^Axin - i 

i=0 



n-l 

I 

i=l 



n-l 

I 

1=1 



This yields vech(a„) = vech(£) +L"=/ v ech(x M _,-^_ ; -) where, for all i 6 N*, ^ is defined by i^,- := 
[B*- l A] , l + [B'- 2 A] j 2 + . . . + [S'^A] j with the convention fi° := /, B' := if i < and [M] t j is 
the <2(<2+ 1)/2 x d(d + 1)/2 block from lines 1 to d(d + 1)/2 and from columns (j - 1 + 1 )/2 + 1 
toj'd(d + l)/2 ofM. 

Thus, W is the Zariski closure of the orbit 

St= (J {?("): x 1 ,...,x„GM rf } 

U | r + y L vBchO^-i^-i))' vech ( x «-p+i-^-/»+i-/)) ? ■ 

V. 



lit t 



n—p 

I 

i=l 



where xi_ 9 = #2-# = ... =xo = (since y(0) = T = (a',0, . . . ,0)')- 

In particular, this implies <p(W n £/ x C W n t/ (cf. (F2)), because <p(£/ x M J ) C U and <p(S r x 
C St yields 



'<p(S T xM d ) C z St= W, 



since <p is a regular map and thus continuous with respect to the Zariski topology. 

Theorem 4.8. Suppose that (H2) holds and that there is a strictly stationary solution (X n ) ne zfor the 
standard GARCH(p,q) model, then the process {Y n )n& takes its values in the algebraic variety of 
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states W fl U. Moreover, one has 

vech(E ;1 ) = (i- £ Bj) vech(C) + £ K t vech(X„_,^_,.) . 
V 7 =i / i=\ 

Proof. Let (X„) ne z be a strictly stationary solution of the standard GARCH(/?,g) model with condi- 
tional covariance matrices E„. We denote by X(n) and E(«) the following random vectors: 

X(n) = (vech(X n X<)V ■ .,vech(X n _ CI+ X-c +i yy 
and E(n) = (vech(Z„) f , . . . , vech(£„_ p+ i) ; )' . 

Since E(n) = ^ +AX(rc - 1) + BE(n - 1) (cf. @~9]>), iterating yields 

/t-i k 
E(n) = £fiVi+B /: r(?i-/t) + £B^ 1 AX(n-/) (4.10) 

;=0 i'=l 

for all icGN. 

Now for any M = (vech(Mi )',..., vech(M p ) f )' and N = (vech(#i )',..., vech(A^)' )' in (R d ( d+i y 2 )P, 
let us denote M > N if and only if Mi >N\,...,M P >N P (where, for all Mi,Nj <E E>d, M,- > N t 
Mj-Ni > O^Mi-Ni positive semi-definite). This defines a partial order on (R rf ( rf+1 )/ 2 )/\ 

Then (l4T0l yields E(n) > J^ 1 fi'^. Since E(n) is finite the series B'^i converges as k -)• oo 
(see for instance [33 1 for further details concerning partially ordered topological spaces; in particular 
the Corollary after Lemma 5 proves that our series must converge). Setting a := TT=o^'^' h * s easv 
to see that a = ^fi + fia. Using the definitions of B and ^ , we obtain that a = (a[ , a{ , . . . , o{) for 
some Oi € 

R d(d+i)/2 which Mfils ffi = vech ( C ) +-£ p j=l BjOi. One may then verify that 0\ = vech(Ei) 

for some Ei € and hence that the spectral radius of YJj=\Bj is l ess than 1 (cf. Proposition I4.3I ). 
Due to Proposition I4.5l (i) we obtain that the spectral radius of B is also less than 1. Thus a = G. 

Next, since the spectral radius of B is less than 1, the sequence (B )k e ^ converges to zero as k — > oo. 
The random vectors (E(n — &))/t e N have a constant law because (X n ) ne z is supposed to be a strictly 
stationary solution of the GARCH model. Thus B k l.(n — k) converges to zero in probability when 
k — > oo. 

With an analog argument as for Y^o^^i one can see that £f =1 fi' _1 AX(n — i) converges almost 
surely as k — >• oo. Hence, taking the limit of (14.101) yields 

oo 

Z(n) = ff + £fl ,-_1 AX(/i-i) a.s. 

i=l 

Using the matrices K{, defined during the investigation of the variety of states W, we obtain 

oo 

vech(E„) = vech(E) + £ Ki vech(X n _ ; X^_ ; ) a.s. 
i=l 

This shows that (Y n ) ne z takes its values in the variety W and hence in W fl U. Note that the strictly 
stationary solution is causal. To finish the proof we refer to Remark 14.61 from which we obtain 
vech(E) = (/ - Ij =1 Bj)- 1 vech(C). □ 
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4.5. Foster - Lyapunov Condition (fFLjl 



We now derive a function V satisfying the Foster-Lyapunov-condition provided that the spectral radius 
of £? =1 Aj + Ylj=\ Bj i s l ess than 1. That is, we prove the following theorem: 

Theorem 4.9. Suppose that the spectral radius of the matrix £f =1 A; + Yfj=\B j * s ^ ess than 1. TTzerc 
?/jere erisf a function V :U — >■ [1,°°) owe? positive constants a < I, b <oo as well as a Borel set K in 
W Dll such that the AFLi - condition is satisfied, i.e. there are positive constants 05 < 1, b < °° such 
that 

PV{x) <a-V{x)+b-l K (x)Vx€WnU. 



Proof. For notational convenience we suppose that in the BEKK representation (12.41 ) /, = Sj = 1 for 
all i = 1, . . . ,q and j = 1,. . . ,p, since the extension to general lj,Sj ^ 1 is obvious and trivial. We set 
A, := A (j i and B ; := S/^. That is, we have 

Zn = C + £A £ X„_^_^ + £ BjZn-jB'j- (4. 1 1) 

i=l 7=1 



If the spectral radius of the matrix Lf =1 A,- + £j =1 fiy i s l ess than 1, then, due to Proposition 14.31 and 
Remarkg31 there exists £ G such that I = C + £? =1 A^£A; + X^ =1 B'j'LBj. 
We define the map V :£/—)• [1,°°) by 

7(F„) := tr^) + . . . +tr(y p £„_ p+1 ) +X% +x X n + ... +X t n _ q+x V p+q X n _ q+x + 1 

where tr( • ) denotes the trace of a matrix and the d x d matrices (V/)i<,-<p+ 9 are given by 

:= £Z^±1 C + f B'j'LBj, l<k<p 
P + q p k 

Setting y = (vech(£„_i) f , . . . , vech(E„_ f7 ) f ,X^_ 1 , . . . ,^_ ? )' G C7 we obtain 

E[y(F„)|F n _! = y] =E [tr^r^+^V^lF^! = y] + tr(y 2 £„_ 1 ) + ...+tr(y p £ n _ p+1 ) 

+x,;_ 1 y p+2 x„_ 1 + . . . +x^ +1 y /7+(? x„_ (?+1 + 1. (4.12) 

Using (14.111 ) for £„, we deduce for the first term at the right hand side 

E [tr(ViZ») +X^ +1 X re |F„_! = y] 
= E [X/y p+1 X„ |F n _! = y] + tr(ViC) + tr^Ai^-i^Ai ) + ... + tr^A,*^-^) 

+ tr(V 1 5 1 £„_ lJ B ; 1 ) + . . .+tr(ViB p L n _ p B t p ) 
= E [X/yp+AlF^j =y] +tx(V 1 C) +X t n _ l A\V l A l X n _ l + ... +X t n _ q A t q V 1 A q X n _ q 

+tr(B\V l B l L n _ 1 ) + . . .+tr(B t p V l B p L n _ p ). 
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Since X n = l}J 2 e n , l}J 2 l}J 2 = L n and E[£„£^] = Id, we obtain 

E [x' n V p +lXn\Yn-i =y] =E[tt(X n (V p+l X n ) t )\Y n _ l =y] 

= tr (E [XXV P+ 1 \Y n -i = y] ) = tr (E [X n X' n \Yn-i = y] V p+l ) = tr(L n V p+l ) 
= tr{V p+l C) + triVp+yA^X-A) + ■■■ + tv{V p+l A q X n - q K-A) 

+ tr(V p+1 B 1 E n _ 1 /? 1 ) + . . . + tx{V p+l B p L n _ p B t p ) 
= tx(V p+1 C) +X t n _ l A\V p+i A l X n _ l + ... +X t n _ q A t q V p+1 A g X n _ g 

+ tv(B\V p+1 B 1 L n _ 1 ) + . . .+tr(B t p V p+1 B p I m _ p ). 

Hence, (14.121) can be rewritten as 

E[y(y„)|y„_ 1 = y ] 

= tr [(B\ (7i + V p+l )B { + V 2 ) E„_i] + . . . + tr [(B^ (Vi + V p+1 + V» E„_ p+1 ] 

+ tr [B^(Vi + V p+ i (A^Vi + V^-ijAi + V p+2 )X n -i + . . . 

+ K- q+ i (4-i( v i+Vi)Vi + %)Wi +^-«A,(Vi + V^i)A^*-, 
+ tx[(V l +V p+1 )C] + l. 

By definition of V}, we deduce 

BU^i + ^ + i)Bt + V fe+ i=V fe --^, 1 <£</?- 1 
B t p (V 1 +V p+l )B p = V p — 



p + q 

A[(Vi + V p+l )A k + V p+k+l = V p+k - , 1 < * < q- 1 
4(V 1+ V, +1 )A, = V, + ,--^-. 

Furthermore, V k is symmetric positive definite for all & = 1 , . . . , p + q which implies that V k — -p^ is 
symmetric positive semi-definite for all k = 1 , . . . , p + q. 
Consider the non-negative constants (a k )i< k < p+q defined by 

a k := max < x 1 [ V k — ) x : x e M. d , x'VkX = 1 i . 

I V J 

Since the maximum is calculated over the unit sphere with respect to the induced norm by V k and 
since this unit sphere is compact, there exists x k £ such that x t k VkXk = 1 and 

t ( C \ t C 

a k =x k [V k — )x k = 1 -4 — — jcj. 

The matrices V k , k = I,. . . ,p + q, and are positive definite which yields < a k < 1 for all k = 
\,...,p + q. 
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Setting Oo := maxfc^ : k = 1 



. , p + q} we obtain < Oo < 1 and 



C 



< OoV k \/k£ {l,...,p + q}. 



p + q 



Hence, for all M G §++ and all k G {1, . 



tr V fc 



C 



) 



M <aotr(V k M). 



We deduce E[V(y„) |7„_i =y}< (XoV(y) + tr(IC) + 1 - Oo. 

If we choose a := (ao + 1)/2 G [1/2, 1) and & := tr(IC) + 1 -a G f0.o°). then the (IFD - condition 
is satisfied with the set K given by 



4.6. Proof of Theorem [231 

Now we can prove our stationarity and ergodicity result for standard GARCH(/?,g) processes. The 
main remaining problem is that K is not compact and, hence, it is somewhat tricky to prove that it is 
small. 

Proof of Theorem \2~4\ . (i) Since due to Proposition 14. 5 1 (ii) the spectral radius of Y?i=\Bj * s a l so l ess 
than 1, Proposition [477] and Theorem 14. 8 1 imply that (A2) and (A3) hold. Using then Proposition 13. 101 
we deduce that (Y n ) ne -^* is i//-irreducible and aperiodic on the state space (W nU,£$(W OU)). 

Define U c = (vech({;e G S j + : x > C})) p x (R d )<? which is a closed set and a proper subset of U. 
Then we have by inspecting the iteration that € W Pi Uc for all Yq^W CiU and & GN,fc>/?.By the 
way, T G Uc by (14.71) . By condition (H3), Theorem 14.91 ensures the existence of a function V which 
fulfils the (IFLI ) - condition on the set K. Now we show that K is small. 

Det K\ = K n U L C and Ki = K\K\ . Using the self-duality of the cone of positive semi-definite 
matrices, it is straightforward to see that V maps unbounded (with respect to norms on (W l ( d+l )l 2 y x 
(M. d ) q ) subsets of U to unbounded subsets of M + and thus K is a bounded subset of U. Inspecting the 
iteration defining the GARCH processes further we see that Y p is not only in W n Uc when Yo G K\, 
but necessarily also in a compact set K C W n Uc conditional on ||e,-|| < f] for i = 1,2, ... ,p and a 
fixed T] > 0. W.l.o.g. one can assume K D K 2 . This implies P p (x,K) > P(\\ei \\ < ri) p =: £ > for all 
leif, due to (H2). 

Moreover, the Markov chain (F M )„ eN * has the Feller property, as an elementary and standard dom- 
inated convergence argument shows, and supp y has non-empty interior (see Proposition I3.101 l. 
Thus, |[26l Proposition 6.2.8] shows that K is petite (see 11251 "). i.e. there is a non-degenerate meas- 
ure v on 8§(W n U) and a probability measure a on N* such that £^Lj a({i})P ! '(jc,B) > v(S) for all 
x £ K and Borel sets B QW C\U . Using Chapman-Kolmogorov this implies 0.5 XILi a({i})P' (x,B) + 
0.5lX= P +\a({i- p})P (x,B) > 0.5Cv(B) for all x eKUKi and Borel sets BQWDU. Thus £u#i 
is petite. Since K C KUK\, also ^ is petite and thus small by GST, Theorem 5.5.7]. 

Applying Theorem l3.11l and Theorem 13. 13 1 we obtain the claimed positive Harris recurrence, geo- 
metric ergodicity as well as geometric /3-mixing and 7i(V) < °° for the stationary distribution %. 




□ 
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Let (X n ) ne z now be the unique stationary GARCH process. Then 7i(V) < °° implies 

E[X' n V p +iXn] < nV(Y n )} = n(V) < oo\/« G Z 

by definition of V (cf. proof of Theorem I4.91 l. This shows that X„ G L 2 for all n € Z. 

Since E[£„] = E[X„X ; '], we deduce E[£„] < °°. Using the diagonal dominance property of a positive 
semi-definite matrix (|mj/| < 0.5(mu +mjj) for M = {mij)\<i,j<d S S^), this implies £„ G L 1 for all 

neZ. 

(ii) We now assume that there is a weakly stationary solution for the standard GARCH (p,q) model. 
Then £ := E[X„X^] is well-defined. Since £ = E[£ n ], taking the expectation in (12.41 ) on both sides 
yields 

1 = C + I I + 11 r £S f ; - r . 

i=U=l j'=lr=l 

Due to Proposition 14.31 the spectral radius of the matrix (J^Lj-A, +Y?j=\ Bj) has to be less than 1. □ 

Note that the proof shows that 71 is concentrated on W n Uc, so in the stationary regime the GARCH 
covariance matrices are always bigger than or equal to C. 

A. Algebraic Geometry 

In this appendix we summarise the necessary details of algebraic geometry to understand the statement 
of our main result. For more details and comprehensive treatments we refer to P]|28l. 

We denote by M.[X\, ... ,X n ] the polynomial ring in n variables formed from the set of polynomials 
in the variables X\, . . . ,X n with coefficients in the field R. 

Definition A. 1. (i) A subset VCR" is called semi-algebraic if it admits some representation of 
the form 

V = {Jf){x£R n : Pij(x) ~ y 0} , 
i=ij=i 

where, for all i = 1 , . . . , s and j = 1 , . . . , r,-, 

(a) ~ y G {>,=,<} 

(b) P iJ (X)eR[X],X = (X l ,...,X n ). 

( ii) A subset l^Ci" is called algebraic if it can be represented as 

V = {x G M" : Pi (x) = . . . = P k (x) = 0} 
where lc£W and P t (X) G R[X X , . . . ,X„}for all i = 1 , . . . ,k. 

Remark A.2. Real algebraic sets can be represented by one single polynomial, namely, if V = 
{/>! = ...= p k = 0}, then we can take P := P\ + . . . +P£. 

Definition A.3 ("Zariski topology"). The topology over R" for which the algebraic sets in R" are the 
closed sets is called the Zariski topology. 

Remark A.4. (i) The Zariski topology is not Hausdorff (i.e. it does not separate points). 
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(ii) Every Zariski closed set in M" is also closed in the usual topology on M". Thus, the usual 
topology is finer than the Zariski topology. 

(iii) We define the Zariski closure of a set A by 7r A := f] B. 

B Zariski closed 
BDA 

Definition A. 5. An algebraic set VCR" is said to be irreducible if it cannot be decomposed as 
V = V\ U V% where both V\ and V2 are algebraic sets and V\^V and V2 ^ V. 
IfV is an irreducible algebraic set, it is also called algebraic variety. 

Definition A. 6. Let V CM" be an algebraic variety and define the ideal of V by 

I(V) := {P € M[X h . . . ,X n ] : P(x) =0 Vx G V} . 

It is an easy consequence of the Hilbert Basis Theorem (cf.for example H23\l ) that the ideal I(V) has to 
be finitely generated, i.e. there exist I S N* and Qi, ■ ■ ■ ,Qi £ I(V) such that I(V) is the ideal generated 
by these polynomials. We then call 

p(V) := suprank ( ^-(x) ) 1 <,</ 

xeV \OXj / l<j< n 

the rank of the ideal I(V). 
A point xq € V is said to be a regular point ofVifp(V) = rank ( ^ (xo) ) 1 < ; < / • Otherwise xq is 

i<j<n 

called a singular point ofV. We write &(V) to denote the set of regular points ofV and J^(y)for the 
set of singular points. 

A natural class of maps are those such that preimages of algebraic sets are again algebraic, i.e. maps 
which are continuous with respect to the Zariski topology. 

Definition A.7. Let V C M" and W C M m be algebraic varieties. Then f :V — > W is said to be a 
regular map, if all its components (f)i<i< m are regular functions, i.e., for all i = 1, . . . ,m, there exist 
Pi, Qi e M[Xi, . . . ,X n ] such that 

P(x) 

Vn{jc€K":gi(*)=O} = and fi(x) = Vx G V. 

Qi \ x ) 



Proposition A.8. Let V C M", W C R m be algebraic varieties and f :V —>W a regular map. Then f 
is continuous with respect to the Zariski topology. 

B. Theory of Markov Chains 

In this appendix we recall the theorems for Markov chains used in the proof of Theorem 13 .121 To this 
end let (X t ) te fq be a Markov chain on the state space (S, &(S)) with transition probability kernel P. 

Theorem B.l (cf. (26J, Theorem 14.2.2). 

Suppose that the non-negative functions V, f, s satisfy the relationship 

PV(x) <V(x)- f(x) + s(x) Vx € 5, 
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then, for each x G S and any stopping time T, we have 



fc=0 



<V(jc) + E, 



T-l 



Jt=0 



Proposition B.2 (cf. E6l . Proposition 9.1.7 (ii)). 

Suppose that (X t ) te ^ is ^-irreducible. If there exists some petite set C G such that L(x,C) = 1 
for all x G 5, then (X t ) te ^ is Harris recurrent. 



Theorem B.3 (cf. Eg, Theorem 14.0.1). 

Suppose that the chain (X t ) te fq is yr- irreducible and aperiodic and let f > 1 be a function on S. Then 
the following conditions are equivalent: 

( i) The chain is positive recurrent with invariant probability measure K and 

n(f) = Jn(dx)f{x) < oo. 

(ii) There exist some petite set C G £§{S), a positive constant b < oo and some extended-valued non- 
negative function V satisfying V(xq) < oo far some xq G S and 

AV(x) :=PV(x)-V(x) < -f(x) + bl c (x)Vx(ES. 



Theorem B.4 (cf. (261, Theorem 15.0.1). 

Suppose that the chain (X t ) te fq is ^-irreducible and aperiodic. If there exist a petite set C G S${S\ 
constants b < oo, j3 > and an extended-valued function V > 1 finite at some xq G S satisfying 

AV(x) < -pv(x) + bl c (x)Vx G 5, 
then there exist r > 1, R < oo such that for any x G {y G S : V(y) < °°} 



^r n \\P"(x r )-7l\\ v <RV(x) 



(B.l) 



n=\ 



with \\}a\\v '■= sup \n(g)\ for any signed measure }l defined on (S,£S{S)). 

g:\g\<V 
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